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Using approximately 1.3 fb~^ of data collected by the DO detector between 2002 and 2006, we 
measure the lifetime of the B^ meson in the — » J/tp/j,^ +XfinaI state. A simultaneous unbinncd 
likelihood fit to the J/'tp + fj, invariant mass and lifetime distributions yields a signal of 881±80 (stat) 
candidates and a lifetime measurement of t{B^) = 0.448lo'o36 (stat) ± 0.032 (syst) ps. 

PACS numbers: 13.20.He, 14.40.Nd, 14.65.Fy 



One of the most interesting mesons that can be studied at the Tevatron is the B^. UnUke most b hadrons, the 
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meson comprises two heavy quarks (& and c) that can 
each decay with significant contribution to the decay rate, 
or they can participate in an annihilation mode. The 
meson is therefore predicted [l|, Q to have a hfetimc of 
only one-third that of the other B mesons, the shortest 
of all weakly-decaying b hadrons. Example final states 
where the c quark acts as a spectator while the h quark 
decays weakly are J/ipn^, B^ J/tpDf, and 

Bf J/ipi^v. 

In this Letter we present a measurement of the lifetime 
of the B^ meson in the B^ J/tpn^ + X final state 
with J/ip A^^M I using approximately 1.3 fb^^ [3] 
of data collected with the DO detector (J] at the Fer- 
milab Tevatron collider. The detector components most 
important to this analysis are the central fiber tracker 
(CFT), the silicon microstrip tracker (SMT), and the 
muon system [Hj. An inclusive muon triggered sample is 
used where events selected only by lifetime-biasing trig- 
gers are excluded. The invariant mass of the resulting 
trimuon system is used to help separate the signal and 
background components and determine their normaliza- 
tions. 

The decay length used to extract the B^ lifetime is 
measured as the distance between the reconstructed pri- 
mary proton-antiproton interaction vertex and the sec- 
ondary vertex formed by the J/ip and the third muon. 
The presence and behavior of the B^ signal is demon- 
strated using mass fits following decay length require- 
ments. We construct models of the lifetime distributions 
of signal and various background components and then 
perform a simultaneous fit to the trimuon invariant mass 
and lifetime distributions to measure the lifetime of the 
B^ meson. 

To simulate B^ properties in this final state and to de- 
termine appropriate selection criteria, Monte Carlo (MC) 
signal samples of B^ J/ip{^ are gen- 

erated using the standard DO simulation chain, includ- 
ing the PYTHIA event generator Q interfaced with the 
EVTGEN decay package followed by full geant [8] 
modeling of the detector response and event reconstruc- 
tion. For the simulated signal samples, the ISGW 
semileptonic decay model 0] for Bf is used. A separate 
sample using a phase-space decay model is generated for 
systematic studies. Another possible decay of the B^ is 
Bf iJj{'^S)n^ + X where '4>{2S) J/V'tt+tt^, and a 
sample of this mode is generated as well. To model one 
of the backgrounds, a large MC sample of inclusive J/ip 
events, including b production via gluon splitting and 
fiavor excitation, is used, only requiring a generator-level 
J/ip ^ M"*"/^" decay. 

We begin by selecting a subsample of events con- 
taining at least one J/ip — > fJ.'^iJ.~ candidate with at 
least two muons of opposite charge reconstructed in the 
CFT, SMT, and the muon system. The track of each 
muon must have transverse momentum pT > 1.5 GeV, 
and match hits in the muon system, or it must have 



Pt > 1.0 GeV and a calorimeter energy deposit con- 
sistent with that of a minimum-ionizing particle. For at 
least one of the muons, hits are required in all three lay- 
ers of the muon detector, and each must have at least 
two hits in the CFT and at least one hit in the SMT. 
The signal region is defined in terms of the dimuon mass 
to be 2.90 < M(^+^^) < 3.26 GeV. The muon momenta 
are adjusted according to a mass-constrained fit to the 
known J/iJj mass [lo| . 

Once a J/V' is found, an additional third track that 
can be associated with the J/tp vertex is sought. The fol- 
lowing cuts are applied to the resulting J/'i/'-l-track can- 
didate: the third track must have at least two hits in 
the SMT, the extrapolation of the three-track momen- 
tum must be consistent with coming from the primary 
vertex, (third track) > 3 GeV, p(third track) > 4 GeV, 
Pt{J/iP + track) > 5 GeV, the probability to form 
a common vertex is greater than 1%, angle between the 
J/t/j and third track < 1 rad, and cos 9 < 0.99 where 9 
is the three-dimensional angle between any two muons. 
If more than one J/4>+tra.ck candidate is present in an 
event, the candidate with the lowest of the J/V'+track 
vertex is selected. To be considered a signal candidate, 
the third track must be identified as a muon: it must have 
hits in all three layers of the muon detector and have tim- 
ing signals in the muon scintillator detectors consistent 
within 10 ns of the beam crossing to reduce contamina- 
tion from cosmic rays. The mass of the J/ ip+n candidate 
is required to be in the range 3 < M{J/ip^) < 10 GeV, 
resulting in a sample containing 14753 events. 

The invariant mass of the J/ip+fi can be used to char- 
acterize and separate each of the components that con- 
tribute to the J/4>+fJ. candidate sample. There are six 
contributions (one signal, and five backgrounds): B^ 
signal (SI); a real J/ip associated with a "fake" muon 
due to a track (JT); fake J/ip mesons from combina- 
torial background (CB); a real muon forming a ver- 
tex with a real J/ip where neither is from a B^ decay 
(JM); B+ J/tPK+ followed by the decay in flight of 
fi'^v; and a cc contribution, where a prompt J/tp 
is associated with a muon (PR). Each component and the 
determination of its mass template is described below. 

The signal mass template is determined from the signal 
MC sample. Theoretical estimates predict the B^ 
J/ipfi^ + X branching fraction to be approximately 5 to 
100 times larger than that of Bf V'(25')/Lt=^ +X [l|,[lH. 
This difference gives 13% to 0% feed-down contribution 
and we take (6.5 ± 6.5)% in the analysis. 

The invariant mass of the J/tp+track in the data sam- 
ple is used to model the JT component. The rate of 
what are denoted fake muons is small and primarily due 
to decays in flight of tt^ fi^v and fi^v. The 

_B+ — > J/ tpK^ decay is used to measure the contribution 
of this component. Fits are made to the i?+ mass peak 
in the J/^p+fi sample and the J/ip+tTa,ck sample, and 
the ratio of the number of B~^ events in the two samples 
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is taken as the fraction of events that are due to a real 
J/^ but a fake ^. Contributions due to B'^ J/'07r+ 
or 5+ J/-!/'7r+ are estimated to be neghgible. 

To describe the CB component, a normahzed mass 
template is formed from events in the J /if) mass side- 
bands. The sideband regions are defined to be events 
with M{^i+n-) in the range 2.62-2.80 GeV or 3.40- 
3.58 GeV. The normalization is taken from the fitted 
number of background events in the signal region under 
the J/ijj mass peak. 

The JM component represents a significant back- 
ground that is dominated by hh backgrounds, where one 
long-lived h hadron decays to J /ip + X and the other de- 
cays semileptonically to a muon (or via a cascade decay 
b ^ c ^ /i). The requirement that the J/ip and /i be 
close in angle increases the relative acceptance for bb pro- 
duction via gluon splitting. To model this background, 
the J/ip QCD MC is used with the requirement that the 
parent of the J/ ip does not arise from a prompt B^ me- 
son, B^ , or cc (the latter two components are estimated 
using the data and described below). 

For the component, a fit is made to the mass peak 
of the B^ in the J/ip+^ data sample. This fitted distri- 
bution is then used as a mass template for the _B+ com- 
ponent, thus reducing the imcertainty in the modeling of 
the width of the mass peak. 

Candidates with L^y < 0, where L^y is the transverse 
decay length defined as the displacement of the J/ip+fi 
vertex from the primary vertex |l2| projected onto the di- 
rection of the transverse momentum vector of the J /^+^ 
system, are used to estimate the mass template of the PR 
component. 

To check the validity of the modeling of the M{J/'ipfi) 
distribution, a fit is first made on the mass distribution 
of the J /ip+^b sample using the templates of the six con- 
tributions described above. Separate fits are made as 
the requirement on L^y is raised to increasingly sup- 
press background. Good agreement of the fitted mass 
components is observed at all Lxy values. To further 
check for the presence of the Bf^ signal, a requirement 
is placed on the transverse decay length significance: 
Lxyl <j(Lxy) > 4, where a[Lxy) is the uncertainty on L^y. 
Figure [U shows the fit to the mass distribution after sub- 
tracting the J /'ip sideband and i?"*" components. In this 
sample, the statistical significance of the Bf signal com- 
ponent is 6.4(7. 

The lifetime of the B^, r, is related to the transverse 

decay length by L^y — ct ■ ^ , where px and m are 

the transverse momentum and rest mass of the B^, re- 
spectively. When the Bf meson decays semileptonically, 
it cannot be fully reconstructed due to the escaping neu- 
trino, and thus pt{B^) cannot be determined. The pT of 
the J /ip + fi system is used instead as an approximation. 
A correction factor, K = pt{J/iPh)/pt{B^), determined 
using signal MC, is introduced to estimate pt{B^). To 
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FIG. 1: Fit to the mass of the J/tp + fi vertex with J/i/i mass 
sideband and B'^ components subtracted and decay length 
significance Lxy/<j{Lxy) > 4 required, (color online) 



obtain the Bf lifetime, the visible-proper decay length 
(VPDL) is measured, defined as L^y-^^^^^^ = ^. 

The if-factor distribution is applied statistically by 
smearing the exponential decay distribution when ex- 
tracting cT{Bf) from the VPDL distribution in the life- 
time fit. The mass of the B^ is taken from To take 
advantage of events with better resolution, the if-factor 
is applied in the analysis in six bins of M{J /ipfx). 

An unbinned likelihood fit is used to measure the av- 
erage lifetime, maximizing C over all i candidates, where 

C = X{[fjTT)r + (1 - .fjT){fcBThB + 

i 

(1 - fcB){fsi^SI + fjM^JM + fB+^B+ + 

{l-fsi-fjM-fB^)rpR})]- (1) 

Each component T consists of a combination of a mass- 
shape template and a lifetime functional model, each de- 
scribed below, to allow for a simultaneous fit of the frac- 
tion components and t{B^). The fractions / of each 
component have been described earlier. The fraction 
fjT = 0.034±0.002 is taken from fits to the B+ peak and 
fcB ~ 0.667 ± 0.004, is found from J/ip mass sideband 
fits. The exponential function J-si is convoluted with a 
Gaussian resolution function and smeared with normal- 
ized X-factor distributions. The width of the Gaussian 
resolution function uses the event-by-event uncertainty 
a{Xi) on the VPDL, multiplied by a floating scale factor 
s to take into account any systematic underestimate of 
cr(Ai) due to tracking systematic uncertainties. A dou- 
ble Gaussian function, centered at VPDL — 0, is used to 
model TpR. The width of the inner Gaussian is given 
by s ■ cr(Ai), and the multiplicative factor for the width 
of the outer Gaussian is determined using MC samples 
and data candidates with negative decay length. Fits are 
made to the respective VPDL distributions to obtain !Fjt 
and J- SB- Empirical functional forms are used for both 
as fixed shapes, and the normalizations via the fractions 
fjT and fsB are allowed to float in the flt. !Fjm consists 
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of a negative-slope exponential and two positive-slope ex- 
ponentials. Starting values for the functional parameters 
are determined from fits to the inclusive J /-ip MC sam- 
ple. The normalization of the negative-slope exponen- 
tial, along with the normalization and slope of one of 
the positive-slope exponentials are constrained by these 
fits. The slope of the negative-slope and second positive- 
slope exponential are allowed to float freely in the final 
fit. is a single exponential function with slope con- 

strained to the world-average value convolved with 
the same Gaussian resolution function as for the signal 
lifetime model. 

Before examining the fit to the data, possible lifetime 
biases are studied. Signal MC samples mixed with back- 
ground are generated with different lifetimes. Fits to 
these samples and ensemble tests indicate no significant 
bias and demonstrate the validity of the extracted statis- 
tical uncertainty. 

A simultaneous fit to the invariant mass and VPDL 
distributions is performed using all the components de- 
scribed above. The fitted lifetime of the meson is 
found to be r(B±) = 0.448lJ5;55i^^ (stat), with an esti- 
mated signal sample of 881 ± 80 (stat) candidates. The 
fitted value of the scale factor is s = 1.35 ±0.02. Figure[2] 
shows the VPDL distribution of the J /ifj+fi sample with 
projections of the fit result overlaid. 
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FIG. 2: VPDL distribution of the J/V'M sample with the pro- 
jected components of the fit overlaid, (color online) 



Stability checks made by dividing the data in half 
based on various selections show no significant lifetime 
variations. The systematic uncertainties considered are 
discussed in detail below and are summarized in Table U 
Variations of the mass within its measurement un- 
certainties [13] make a negligible difference in the life- 
time. The B^ signal modeling uncertainty is estimated 
from the difference between the default and phase space 
decay models. The uncertainty in pt{B^) is found by 
reweighting the spectrum to correspond to varying the 
factorization and renormalization scales iip =m,R = 



-\/py(parton) + by factors of a half and two lj| . To 



address uncertainties on the predicted prib) for the sig- 
nal and background component distributions, a momen- 
tum weighting factor that is applied to MC samples to 
improve the simulation and include the effects of the trig- 
gers is removed. In all of the above cases, both new signal 
mass templates and i^-factor distributions are generated 
and the analysis repeated. To assess the systematic un- 
certainty due to the modeling of the inclusive J/tp MC 
mass distribution, contributions due to bb production via 
gluon splitting and then flavor excitation are entirely re- 
moved. The shape of the mass template for the prompt 
component is varied within the statistical errors of its de- 
termination, and the observed lifetime variation assigned 
as a systematic error. All of the systematic uncertainties 
described above are added in quadrature and summarized 
in Table 1 under the category of mass model uncertain- 
tics. 

To test the assumption that the modeling of the life- 
time of the J/ip combinatoric background can be approx- 
imated by taking the average of the upper and lower mass 
sidebands, the fit is performed using only the high or 
the low mass sideband, and a systematic uncertainty of 
one half the resulting shifts in lifetime is assigned. The 
scale factor s is varied over the range of values, 1.2 — 1.4, 
observed in other lifetime analyses [l3| as well as as- 
signed a functional form and the variation in lifetime 
assigned as a systematic uncertainty. In the modeling 
of the prompt lifetime PDF, a single Gaussian function 
rather than a double Gaussian is used to describe the zero 
lifetime events. The shape parameters of the sideband 
lifetime model are changed by varying the fit parame- 
ters within their uncertainties. The parameters defining 
the J/tp QCD MC lifetime model are varied around their 
central values by ±lcr. For the B'^ lifetime model, the 
central value is changed by zLla. The B^ lifetime is also 
allowed to fioat as a systematic study on the B^ lifetime 
as well as a check of the JS"*" lifetime, finding a value of 
1.88 ± 0.19 ps, consistent with the world average value 
of 1.638 ± 0.011 ps [lO|. All of the systematic uncertain- 
ties described above are added in quadrature under the 
category of lifetime model uncertainties. 

Smaller systematic uncertainties arise from the vari- 
ation of the fraction of the feed-down Be ip{2S)X 
signal component between 0% and 13%, and from possi- 
ble alignment effects, estimated using signal MC with a 
modified detector geometry within the alignment toler- 
ances. 

In summary, using approximately 1.3 fb^^ of data, 
the lifetime of the Bf meson is measured in the B^ 
J/tpfi^ -I- V final states. Using an unbinned likelihood si- 
multaneous fit to the J/ijj+ii invariant mass and lifetime 
distributions we measure 

T{Bt) = 0.448l^;^i^^ (stat) ± 0.032 (sys) ps. 

This measurement is consistent with theoretical predic- 
tions of 0.55 ± 0.15 ps in an operator product expansion 
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TABLE I: Summary of estimated systematic uncertainties. 



Systematic source 


Ar (ps) 


Mass model uncertainty 


±0.021 


Lifetime model uncertainty 


±0.022 


Signal feed-down fraction 


±0.005 


Alignment 


±0.006 


Total 


±0.032 



calculation T| and 0.48±0.05 ps using QCD sum rules [li- 
lt is also consistent with the most recent measurement 
from the CDF collaboration [l^, but with significantly 
better precision, making this measurement of t{B^) the 
most precise to date. 
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